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Abstract 

A novel continuum theory of two-dimensional quantum gravity, based on a ver- 
sion of Causal Dynamical Triangulations which incorporates topology change, 
has recently been formulated as a genuine string field theory in zero-dimensional 
target space (arXiv:0802.0719). Here we show that the Dyson-Schwinger equa- 
tions of this string field theory are reproduced by a cubic matrix model. This 
matrix model also appears in the so-called Dijkgraaf-Vafa correspondence if the 
superpotential there is required to be renormalizable. In the spirit of this model, 
as well as the original large- A expansion by 't Hooft, we need no special double- 
scaling limit involving a fine tuning of coupling constants to obtain the continuum 
quantum-gravitational theory. Our result also implies a matrix model represen- 
tation of the original, strictly causal quantum gravity model. 
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1 Introduction 



Dynamical triangulations (DT) were introduced as a regularization of the Polyakov 
bosonic string and of two-dimensional quantum gravity [1, 2, 3]. Using this regu- 
larization, one could show that a tachyon-free version of Polyakov's bosonic string 
theory does not exist in target space dimensions c? > 1 [4]. However, when viewed 
as a theory of 2d quantum gravity coupled to matter with central charge c < 1, 
the theory (non-critical string theory) did make sense. Using matrix-model tech- 
niques and other combinatorial methods, it was sometimes even advantageous to 
use the regularized theory for analytic calculations. Related attempts to use DT 
as a regularization of higher- dimensional quantum gravity [5] were less successful 
[6]. This triggered the introduction of Causal Dynamical Triangulations (CDT), 
which use causal, Lorentzian instead of Euclidean curved spacetimes as a fun- 
damental input. Evidence has been accumulating that they provide us with a 
non-trivial theory of quantum gravity in four dimensions [7, 8]. 

While the higher-dimensional DT and CDT theories of quantum gravity at 
this point rely strongly on numerical simulations, the 2d CDT theory of quantum 
gravity can be solved analytically [9], like its 2d Euclidean DT counterpart. This 
is described in detail in two recent papers, where we have also developed a com- 
plete string field theory in a zero-dimensional target space for the CDT version 
of 2d quantum gravity [10, 11].^ This string field theory or third quantization 
of 2d quantum gravity uses the formalism already developed by Ishibashi, Kawai 
and collaborators for the DT version of 2d quantum gravity in the context of 
non-critical string theory [12, 13]. For non-critical string theory, it is known from 
[12] that the string field theory reproduces the results of the double-scaling limit 
of the matrix models whenever the results can be compared. 

Given the formal similarity between the CDT string field theory and the non- 
critical string field theory, it is natural to ask whether there also exists a matrix 
model which reproduces the results of the former. Below we will show that the 
answer is in the affirmative. However, since the scaling found in the CDT model is 
different from the conventional double-scaling limit of matrix models, a different 
limit needs to be taken. We will show that the limit is simply the conventional 
limit used in the context of the Dijkgraaf-Vafa duality to U (N) supersymmetric 
gauge theories [14]. 

^Due to the inclusion of higher-genus surfaces, this amounts to a non-trivial generalization 
of the original, strictly causal CDT formulation. 
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2 CDT string field theory 



We have recently developed a string field theory for Causal Dynamical Triangu- 
lations [11]. The starting point of the CDT quantization of gravity is the as- 
sumption that in a gravitational path integral over spacetimes with a Lorcntzian 
signature only causal geometries should be included, an idea dating back at least 
to [15]. How this can be done explicitly in a regularized theory, how one can 
rotate to Euclidean signature to perform explicit calculations, and eventually 
take the cut-off (or lattice spacing) to zero is described in detail in [9] for two 
and in [16] for three and four spacctime dimensions. We demonstrated in [10] 
how one can still solve the 2d model analytically when the original formulation 
is extended to allow the light-cone structure to become degenerate in isolated 
points. In [11] we generalized these results to a genuine string field theory, which 
enabled us in principle to calculate the amplitudes of certain spatial correlators, 
for two-dimensional worldsheets of any topology^. 

Let us briefly deflne the CDT string fleld theory, while referring to [11] for 
details. We will work in a Euclidean notation, which means that we started 
out with a Lorentzian signature, regularized the theory, rotated it to Euclidean 
signature as described in [9] and then took the lattice cut-off a to zero. In 
particular, this implies that all quantities discussed below are already continuum 
quantities. 

We have a "free" Hamiltonian Hq which describes the causal propagation of 
a spatial universe with respect to proper time t. Let a spatial universe with the 
topology of a circle of length I2 (the "exit" loop) be separated a geodesic distance 
t from another spatial loop of length l-i (the "entrance" loop), and denote the 
corresponding amplitude by G^^\li, l2',t). It is represented by the path integral 

G^;\hMt)^ I P[5^.] e-^[^-l, (1) 

with the (Euclidean) gravity action 

S[9,..] = A y d'C^det 5^,(0 +xjdh + yj d/2, (2) 

where A is the cosmological constant, x and y are two so-called boundary cos- 
mological constants, g^^ is a metric representing the geometry (diffeomorphism 
equivalence class) [g^u] , which is assumed to be strictly causal in the sense of the 
original CDT model [9]. This means essentially that the topology of its spatial 
sections will not change as time advances. We choose the spacetime to have the 

^For earlier results in this direction we refer to [17]. 
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topology of a cylinder, x [0, 1] . In a Hilbert space language one has [9] 



G^^\h,l2;t) = (/2|e-*^°«|/i), Ho{l) = -i^ + ^i- (3) 

Next, we will generalize the class of geometries integrated over in the path integral 
(1). As a function of time t, spatial geometries will be allowed to branch into 
disconnected circles, and the resulting baby universes can subsequently merge 
again. Furthermore, a spatial universe will be allowed to vanish into the "vacuum" 
if it has length zero. These topology-changing processes can be described by the 
string field Hamiltonian 

H ^ jj ^if\l)Homil) - g j dh j dl2^\h)¥{l2)^{h + l2) (4) 

-ag I dhj dl2^\h + l2)^{l2)^{h)- j jSmil). 

The operator ^ is a "second quantized" Hamiltonian in the sense of many-body 
theory. We introduce creation and annihilation operators ^^{l) and ^{l) for 
universes of length I which act on the above-mentioned vacuum state |0), with 
defining relations 

10 = vi>t(/)|o), vi/(/)|/) = |o), vi/(/)|o) = (oi^t(;) = 0, = m - 1')- 

(5) 

In (4), is a coupling constant of mass dimension 3, and a is a dimcnsionless pa- 
rameter allowing us to distinguish between the merging and splitting of universes, 
which will be set to 1 at the end of the calculation. For a = 1, is hcrmitian 
except for the presence of the tadpole term proportional to 5{l). It tells us that a 
universe can vanish when it has zero length, but cannot be created from nothing. 
Also the two interaction terms have a straightforward geometric interpretation. 
The first term replaces a single spatial universe of length li + 12 with two spatial 
universes of length li and I2, while the second term represents the time- reversed 
process where two spatial universes merge into one, again without changing the 
total length li-\-l2. The couphng constant g clearly takes on the role of string 
coupling constant, since the splitting of spatial universes is associated with a 
factor g and the merging with a factor ag^ making for a combined factor of ag^ 
whenever the spacetime topology changes (sec [11] for a detailed discussion). 

We can use the string field theory associated with H to calculate connected 
multi-loop correlators defined by 

w{h, ...,ln)^ lim (0| e-*^*t(^^) . . . ^t(;^) |0),„,,rf. (6) 

c— >oo 

They describe all possible ways in which an initial state of n spatial loops can 
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Figure 1: A typical geometry in the string field theory contributing to the amplitude 
w(li, . . . ,ln) of eq. (6). Proper time progresses upwards. The dots mark singular points 
of the causal structure. 



evolve and eventually vanish into the vacuum, while forming a connected two- 
dimensional geometry (c.f. Fig. 1). The amplitudes w{li, . . . ,ln) are determined 
from the string field theory "partition function" 

Z(J) = lim (0| e-*^ei''^'-^W*'»|0) (7) 



through the prescription 

d"F(J 



will, ...Jn) 



SJ{h)---SJ{Q 



F{J) = log Z{J). (8) 



In [11] we derived the Dyson- Schwinger equations for the correlators . . . , 
They follow from the ^-independence of Z(J), which leads to the relation 



6^F{J) 



5J{l')5J{l - V) 

One obtains the Dyson-Schwinger equations for the amplitudes w{li, . . . ,ln) by 
differentiating (9) n times with respect to J{1) and then setting J(/) = 0. The 
general equation at order n can be written down easily, but is involved. We will 
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give only the first three equations exphcitly, from which the general structure 
should be clear. The Dyson-Schwinger equations are most conveniently formu- 
lated in terms of the Laplace-transformed amplitudes 

1 /'OO POO 

w{x,,...,Xn) = —^ / dh--- dlne-^'''--^-'-w{h,...,ln), (10) 

a" Jo Jo 

where for convenience we have rescaled the amplitudes by a factor a"'"^, compared 
with the convention used in [11]. Introducing the notation 

V'{x) = ^{X-x^), V{x) = ^(Xx-^x^^, (11) 

we obtain from (9) (see [11] for details) the equations 

= dx(^-V'{x)w{x) +w^{x) + aw{x,x)^ (12) 

= dx(j—V'{x)+2w{x)]w{x,y) + aw{x,x,y)^ + 

dy {[-V'{y) + 2w{y)]w{x, y) + aw{x, y, y)) + (13) 

V X — y / 

= dx(j—V'{x) + 2w{x)]w{x,y,z) + aw{x,x,y,z)^ + 
9y {[-^'(y) + '2w{y)]w{x, y, z) + aw{x, y, y, z)^ + 

9z {[-V'iz) + 2w{z)]w{x, y, z) + aw{x, y, z, z)) + (14) 
2dx[w{x, y)w{x, z)] + 2dy[w{x, y)w{y, z)] + 2d^[w{x, z)w{y, z)] + 
2 /^gi^g^ ^(^'^)~'^(^'^) I I g^gT^^'^^"'^^^'^'' 



X — y X — z y — z 

Let us introduce the expansion^ 



WiXi 



, • • • , a^n) = a^Wh{xi, ...,Xn). (15) 



h=0 



As shown in [11], h can be interpreted as the number of handles of the world- 
sheet, and the equations above can be solved iteratively in h. More precisely, the 



^Note that both w and Wh are still ^-dependent, although we do not write the dependence 
explicitly here. 
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equations at order q;° allow us to determine Wo{x), Wo{x,y), and similarly the 
equations at general order determine Wh{x), Wh{x,y), etc. For example, one 
finds 

M^) = I {V'i^) + -(^ - c)V(a;-c_)(a;-c+)), (16) 



. ^ 1 1 / xy-^(c_ + c+)(x + y) + c_c+ . 

wo{x, y) = --, . ^ — = - 1 , 17 

2 {x - yf \y/{x - c_)(x - c+) V(y - c_)(y - c+) 

where the constants c, c± are determined by 



c 



,3 



Ac + ^ = 0, c± = -c±^/Vc. (18) 



Writing the amplitudes in this fashion leads one to the surprising realization that 
wq{x) and wq{x, y) coincide with the large- iV limit of the resolvent and the planar 
loop-loop correlator [18, 19, 20] of the Hermitian matrix model with potential 

V{M) = ^M - ^M^\ (19) 

This is a potentially exciting result, because so far no standard formulation in 
terms of matrix models has been found for a CDT model, in contrast to the "old" 
Euclidean DT models. We will in the following section prove a more general result, 
which will identify the Dyson-Schwinger equations derived above with the loop 
equations of a Hermitian matrix model with the cubic potential (19). 

3 Matrix loop equations 

Let M denote an N x N Hermitian matrix, V{M) a potential of the form 

oo 

\/(M) = -v^M^ (20) 

k=l ^ 

and define the functions 

W(.., ...,.,) = ^(t. -1^) . . . (tr -i^))^ . (21) 

The subscript c in (Oi(M) ■ ■ ■ On{M))c denotes the connected part of the expec- 
tation value, which itself is defined as 

iOAM) . . .0„(M)) = J 'yj^ . (22) 
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It is well known that the matrix integrals corresponding to (21) possess a large- 
expansion. Assume we have the so-called one-cut solution related to this 
expansion. The invariance of the matrix integral under a change in variables 
leads to the loop equation [18, 19, 20] 

where the integration contour C encloses the cut, but not the point x. From this 
equation one can obtain the equations for the multi-loop correlators by differen- 
tiating with respect to the coupling constants gk in terms of the so-called loop 
insertion operator [19, 21] according to 

'^'-'■■■■■-")- dKfe)''.'."'dKM 

where the insertion operator is given by 

^ ^ (25) 



dV{x) ^x'^+id^fc- 



For a given potential with fixed coupling constants one uses these relations in 
the following way. Assume that gk can vary, act with the loop insertion operator 
on (23) as many times as needed, and then set gk — gt- This leads to the desired 
loop equations. In order to compare with the Dyson-Schwinger equations of our 
string field theory, we differentiate the equations obtained with respect to x, and 
finally find for the potential (19) the equations 



- d:,(^-V\x)W{x) + W\x) +—W{x,x))--, (26) 
= d,([-V'{x) + 2W{x)\W{x,y) + ^W{x,x,y)) + 

\ X — y / 

= d,([-V\x) + 2W{x)]W{x,y,z) + ^W{x,x,y,z)) + 

2d,{w{x,z)W{x,y))+ (28) 

^W{x,z)-W{y,z)^ ^ ^^^^ rW{x,y) -W{z,y)^ 
V X — y / V X — z / 

Using that W{xi, . . . ,a;„) is a symmetric function of its arguments, we see that 
eqs. (26)-(28) lead to exactly the same coupled equations for W as do (12)-(14) 
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for w if we identify 

« = (29) 

In this case the discussion surrounding the expansion (15) is nothing but the 
standard discussion of the large- expansion 

oo ^ 

W{xi, ...,Xn)^Yl 'mh ^h{xi, ...,Xn) (30) 



h=0 



of the multi-loop correlators (see, for instance, [21] or the more recent papers 
[22, 23]). The iterative solution of these so-called loop equations is uniquely 
determined by Wq{x) (and the assumption that W{xi, . . . , Xn) is analytic in those 
Xi that do not belong to the cut of the matrix model) , and we have already seen 
that Wq{x) — Wq{x). 



4 Discussion and Outlook 

Let us consider the matrix model corresponding to the potential (19). We can 
perform a simple change of variables M — > —M — \/A in the matrix integral to 
obtain a standard matrix integral 

Z{m,g)^ J dMe-^^(^\ (31) 

where the new potential (up to an irrelevant constant term) is given by 

V{M) = - + ^M^) , m = 2\f\. (32) 

It is amusing to note that the matrix integral (31) is precisely the kind of matrix 
integral considered in the so-called Dijkgraaf-Vafa correspondence [14], where 
V{(^) is the trcc-lcvcl superpotential of the adjoint chiral field $, which breaks 
the supersymmetry of the unitary gauge theory from J\f — 2 to J\f — 1. If 
one demands that this tree-level potential correspond to a renormalizable theory, 
its form is essentially unique, and precisely of the form (32) originally used by 
Dijkgraaf and Vafa, with g a dimension-three coupling constant coming from 
topological string theory and in the DV-correspondence related to the glueball 
superfield condensate in the gauge theory. 

In the "old" matrix model representation of non-critical strings and 2d grav- 
ity one had to perform a fine-tuning of the coupling constants in order to obtain 
a continuum string or quantum gravity theory. This implemented the gluing of 
triangles (or, more generally, of squares, pentagons, etc.) which served as a reg- 
ularization of the worldsheet. The fine-tuning of the coupling constants reflected 
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the fact that the hnk length of the triangles (the lattice spacing of the dynamical 
lattice) was taken to zero in the continuum limit. The situation here is different. 
Although CDT can be constructively defined as the continuum limit of a dynam- 
ical lattice, we have in the present work been dealing only with the associated 
continuum theory. Thus in our case the matrix model with the potential (19) (or 
(31)) already describes a continuum theory of 2d quantum gravity. Its coupling 
constants can be viewed as continuum coupling constants and the role of is 
exactly as in the original context of QCD, namely, to reorganize the expansion in 
the coupling constant g. 't Hooft's large- A" expansion of QCD is a reorganization 
of the perturbative series in the Yang-Mills coupling gyM, with 1/A^ taking the 
role of a new expansion parameter. In this framework, after the coefficient of the 
term l/N"^^ of some observable has been calculated as function of the 't Hooft 
coupling gjj = gyu^^ one must take A" = 3 for SU (3), say. The situation in CDT 
string field theory is entirely analogous: starting from a perturbative expansion in 
the "string coupling constant" g (in fact, in the dimensionless coupling constant 
g/X^/^, as described in [10, 11]), we can reorganize it as a topological expansion 
in the genus of the worldsheet by introducing the expansion parameter a. For 
the multi-loop correlators this expansion is exactly the large- expansion of the 
matrix model (19) and the coefficients, the functions Wh{xi, . . . , a;„), are exactly 
the multi-loop correlators for gemis-h worldsheets of the CDT string field theory 
with a = 1. 

As a "bonus" for our treatment of generalized (and therefore slightly causality- 
violating) geometries, we also obtain a matrix formulation of the original two- 
dimensional CDT model proposed in [9], where the spatial universe was not al- 
lowed to split. Working out the limit as 51 —> of the various expressions derived 
above, we see that this model corresponds to the large- limit of the matrix 
model where the coupling constants go to infinity, but at the same time the cut 
shrinks to a point in such a way that the resolvent (or disk amplitude) survives, 
that is, 

wo{x) -= = wcdt{x). (33) 

X -\- V A 

The existence of a matrix model describing the algebraic structure of the Dyson- 
Schwinger equations leads automatically to the existence of Virasoro-like opera- 
tors n> —1 [18, 19], which can be related to redefinitions of the time variable 
t in the string field theory. This line of reasoning has already been pursued by 
Ishibashi, Kawai and collaborators in the context of non-critical string field the- 
ory. It would be interesting to perform the same analysis in the CDT model and 
show that reparametrization under the change of time-variable will reappear in 
a natural way in the model via the operators L„. The results should be simpler 
and more transparent than the corresponding results in non-critical string field 
theory since we have a non-trivial free Hamiltonian Hq in the CDT model. 
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